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Abstract 

We examine a simple model of proton pumping through the inner membrane of mitochondria in 
the living cell. We demonstrate that the pumping process can be described using approaches of 
condensed matter physics. In the framework of this model, we show that the resonant Forster-type 
energy exchange due to electron-proton Coulomb interaction can provide an unidirectional flow 
of protons against an electrochemical proton gradient, thereby accomplishing proton pumping. 
The dependence of this effect on temperature as well as electron and proton voltage build-ups are 
obtained taking into account electrostatic forces and noise in the environment. We find that the 
proton pump works with maximum efficiency in the range of temperatures and transmembrane 
electrochemical potentials which correspond to the parameters of living cells. 

PACS numbers: 87. 16. Ac, 87.16 Uv, 73.63.-b 
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I. INTRODUCTION 



A living cell can be considered as a tiny electrical battery with a transmembrane potential 
difference of order —70 mV (with a negatively charged interior). Even a higher potential, 



AV ~ —200 mV, is applied to the inner membrane o: 



a mitochondrion, an organelle, which 



. naa. 



To create and maintain such 



produces most of the energy consumed by the cell, 
an electrical potential, mitochondria employ numerous proton pumps converting energy of 
electrons into an electrochemical proton gradient that is harnessed thereafter to drive the 
synthesis of adenosine triphosphate (ATP) molecules. Translocation of protons across the 
inner membrane of mitochondria is performed by the enzyme cytochrome c oxidase (COX). 
Although crystal structure of COX is known in detail, a molecular mechanism of the redox- 
driven proton pumping remains a mystery despite of the significant latest advances based 
on time-resolved optical and electrometric measurements j^, 5]. 

The electron transport chain of COX consists of four metal redox centers, Cu^, heme a, 



heme as, and Cub 



fly, El. 



The process starts when the mobile electron carrier, cytochrome 
c, moving from the positively charged P-side of the membrane, donates a high-energy electron 
to a dinuclear copper site, Cu^ (see Fig.l). After that, the electron proceeds to the heme a 
with a subsequent transfer to the binuclear center formed by heme and a copper ion Cub, 
where the dioxygen molecule O2 is reduced to water. To produce two molecules of water in 
the catalytic cycle with four electrons (e~) 

O2 + 8 H+ + 4 e" ^ 2 H2O + 4 H+, 

the cytochrome oxidase consumes 4 substrate (chemical) protons which are translocated 
from the negative N-side of the inner mitochondrion membrane to the binuclear center. In 
the process, four more protons (B.^) are taken from the N-side and pumped to the positive 
side (Hp). Here, subscripts N and P for the protons denote the location of the proton H+ at 
the negative (N) or positive (P) side of the membrane, respectively. A residue E278 (for the 
Paracoccus denitrificans enzyme) or a conserved glutamic acid, Glu2A2 (for the bovine 
enzyme , located at the end of the so-called D-pathway llOy, can serve as starting 

points for both substrate and pumped protons on their way from the N-side to the binuclear 
center. In the next phase, a proton is transferred to an unknown yet protonable pump site 
X which is located on the P-side of the heme groups and electrostatically coupled to heme a 
and to the binuclear iron-copper center os/Cub fl S]. On the final stage, the proton moves 



from the site X to the positive side of the membrane after uphill pumping. In the context 
of a pure electrostatic model proposed in Refs. [4, [Sj, the protonation of the site X leads to 
the equalization of electron energy levels in hemes a and that facilitates a transfer of an 
electron from heme a to the binuclear center. This electron attracts a substrate proton which 
moves from the N-side of the membrane to the site X, expelling the first, pre-pumped proton 
to the P-side. Detailed density functional and electrostatics studies of this and other models 
have been performed in sj. 
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15| . However, a mechanism of energy transmission 



from electrons to protons resulting in an unidirectional translocation of protons against the 
concentration gradient is still uncertain. For better understanding of this phenomenon, it is 
useful to combine a comprehensive analysis of the energetic and spatial structure of enzymes 
with simple and physically transparent models. 

In the present paper, we approach the problem taking into account the similarity of the 
electron-driven proton transfer to the quantum transport of electrons through nanostructures 

n 

|16l |. The interaction between electrons and protons is described by a Coulomb potential, 
but, in addition to the standard electrostatic terms, we analyze effects of the Forster-type 
Coulomb exchange on the resonant energy transduction between electron and proton 
subsystems. Each of the subsystems is supposed to have two active sites: Ig, 2e for electrons, 
and Ip, 2p for protons. We consider here the possibility when both electron sites belong to 
the same potential well, localized in the binuclear center os/Cub, while both active proton 
states 2p and Ip can be ascribed to the pump center X (see Fig.l). This positioning of active 
sites corresponds in some sense to the electrostatic model of Ref. jsl, based on time-resolved 
measurements of electron transfer in COX enzyme 4|. 

During the Forster process, an electron moves from the state 2e, which has a higher 
energy, to the state Ig, with a lower energy; whereas a proton jumps from the lower-energy 
state Ip to the higher-energy state 2p (see Fig. 1). The same mechanism is responsible for 
the Fluorescence Resonant Energy Transfer (FRET) in biological systems 13], as well as for 
the exciton transfer in condensed matter 19 1. 

The Forster term originates from the matrix element of the Coulomb electron-proton 
potential between the overlapping wave functions of the electron states 2e and le, and the 
overlapping wave functions of proton states Ip and 2p [20] . Calculations show that this term 
is directly proportional to the product of the dipole moments of electron and proton two-level 
systems, also inversely proportional to the cube of the distance between the electron and 
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proton sites, and requires to satisfy resonant conditions for the energies of the electron and 
proton subsystems. Accordingly, the Forster term is much weaker than standard electrostatic 
terms. However, as a consequence of its overlapping origin, this term opens a new channel 
for simultaneous tunneling of electrons and protons, in addition to the direct tunneling. We 
demonstrate that it is the Forster-type coupling that results in an effective electron-proton 
energy transfer, followed by the proton pumping from the negative to the positive side of 
the inner mitochondria membrane. 

The rest of the paper is structured as follows. Formulation of Hamiltonians and energetic 
spectra of the problem is presented in Section II. Expressions for electron and proton currents 
arc obtained in Section 111. In Section IV, we derive equations of motion for the density 
matrix. In Section V, these equations are solved numerically and the obtained dependencies 
of the proton current on temperature, electron and proton voltage build-ups, and deviation 
from the resonant conditions are discussed. Section VI contains our conclusions. 



II. MODEL FORMULATION 

Electrons and protons on sites u = 1, 2 are characterized by the Fermi operators a+, Oo-, 
and b^,b„, respectively, with the corresponding populations, n„ = a^a„ and N„ = b^b„ (we 
interchangeably use the notation "site" = "state"). We assume that each electron site or 
proton site can be occupied by a single particle, so the maximal populations can be, at most, 
one electron on each one of the two separate electron sites, and, at most, one proton on each 
one of the two separate proton sites. To describe the continuous flow of carriers through the 
system, we assume that the electron site 2 is coupled to the left (L) reservoir, which serves 
as a source of electrons, and the electron site 1 is coupled to the right reservoir (R) playing 
the role of drain. At the same time, the proton site 1 can be populated when protons jump 
from the reservoir located on the negative (N) side of the membrane. On the positive side 
of the membrane, there is another proton reservoir which serves to depopulate of the proton 
site 2 (see Fig. lb). In the framework of this model, here we neglect the couphngs between 
the electron site 1 and the reservoir L, and between the site 2 and the reservoir R. We also 
neglect the tunneling between the proton site 1 and the positive side of the membrane (P), 
as well as the tunneling between the proton site 2 and the negative side of the membrane 
(N). 
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The electrons in the reservoir (lead) a (a = L, R) or the protons in the reservoir (lead) 
f3 {(3 = N, P) can be characterized by additional parameters k and q, respectively, which 
have meanings of wave vectors in condensed matter physics. To describe the electronic and 
protonic sources and drains, we introduce the electron creation and annihilation operators in 
the a-lead as c^^, Cka, and their proton counterparts for the /?-lead as c/^, (ig/3. The number 
of electrons in the a-lead is determined by the operator J2k ''^ka, with Uka = C^a^ka, whereas 
the proton population of the /3-lead is given by the operator J2q with A^^^ = d^pdqjj. It 
is well-known that in real biological structures, couplings between the active sites 1,2 and 
the reservoirs can be mediated by many bridge states, similar to the CuA-site and heme a, 
which can be subjected to conformational changes. Conformation changes can also provide 
a selectivity in coupling between the active sites and the leads [2I. 



A. Electron and proton Hamiltonians 

The Hamiltonian of the electron-proton system incorporates a term related to eigenener- 
gies ea \ of electrons and protons, respectively, located on the sites a = 1,2, as well as 
a term describing electron and proton energies eka, Ekp of the leads a = L, R; P = N, P : 

o" ka qf) 

The Hamiltonian -ffdir, 

i/dir = -Ka^ai - Alata2 - A.bth - Albib2, (2) 

is responsible for the direct tunneling of electrons and protons between the corresponding 
sites 1 and 2, with the rates and A5. Notice that the direct tunneling has a highly non- 
resonant character since the energy levels of the sites 1 and 2 are well separated: e^2^ — ef'^ 3> 
Aa, eI^^ ~E^ ^ A5. To take into consideration the coupling of the active sites 1 and 2 to the 
corresponding reservoirs of electrons and protons, we introduce the tunneling Hamiltonian 

-f^tun = - ^ tkRctjitti - ^ tkL4j^a2 - ^ Tg7vdjy6i - ^ Tqpd'^pb2 + h.c. (3) 

k k q q 

The Coulomb force plays the most important role in the process of energy transfer from 
the electron subsystem to protons. This interaction is determined by the Coulomb potential 



where Ve^Vp are the electron and proton positions in their local frame of reference, and R 
is the distance between the electron and proton sites, R ^ rg, Vp. A direct electron-proton 
Coulomb attraction is determined by the energies u^a-' (c = le, 2e; <j' = Ip, 2p). In addition, 
we take into account the repulsion of the two electrons located at the sites Ig and 2e (energy 
scale ~ Me) jointly with the repulsion of two protons localized on the sites Ip and 2p (an 
energy parameter Up). It should be noted that all energy characteristics Uaa',Ue,Up are 
modified compared to their original values because of Coulomb interactions between the 
active sites and the electron and proton reservoirs. As a result, the Hamiltonian related to 
the direct Coulomb interaction has the form 

= - ^ u^^m^N^, + Uenin2 + UpNiN2. (5) 

crcr' 

B. Forster term 

The direct Coulomb coupling between electrons and protons should be complemented by 
the Forster term. 

Hp = Vpa^a26^6i + V^a^ai6^625 (6) 
which originates from the cross matrix element of the Coulomb potential (jlj) 

This matrix element is taken over the electron-proton wave function | le2p) , with the electron 
being in the state le and the proton being in the state 2p, and the wave function \2elp), with 
the electron being in the state 2e and the proton being in the state Ip. The Forster term 
can be significant in the case of an electron-proton resonance when the distance between 
the electron energy levels ei and €2 is close to the separation of the proton energy levels 
El and E2 : €2 — — E2 — Ei. Therefore, the states |le2p) and |2elp) have almost the 
same energy: ei + E2 ~ e2 + -Ei, that is favorable to transitions between these states. The 
contributions of the other cross-elements of the electron-proton Coulomb attraction, such as 
(2elp|u(re, Vp, _R)|le2p), (2e2p|u(re, Tp, R)\le2p) , ctc, which have a non-resonant character, are 
quite small (~ Vp / {E2—E1) ^ 1 at E2—E1 ~ 500 meV, Vp ~ 1 meV), and can be neglected. 
We consider here a situation where the wave functions le, 2e represent the ground and the 
first excited state of the electron in a parabolic potential well which is placed a distance 
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R from the proton potential well containing two proton states Ip, 2p. Using the expansion 



|r| < |R|), 



R 



1 
R 



i?2 



(8) 



we find that the matrix element Vp characterizing the strength of the Forster term is pro- 
portional to the product of the dipole moments, ero and cRq, of the electron and proton 
sites 1 and 2 and inversely proportional to the cubic power of the distance R between these 
sites: 

= (9) 



27ieoer R^ ' 

For a protein with a dielectric constant = 3 and the electron/proton wave function 
spreadings tq = 0.1 nm and Rq = 0.01 nm, we estimate the Forster matrix element as 
Vf ~ 1 meV, if the distance between the electron and proton sites R = 1 nm. 



C. Dissipative environment 



To account for the effects of a dissipative environment on the electron and proton transfer, 



we resort to the well-known model 



23 



2J] where the polar medium surrounding the 



electron and proton active sites is represented by two systems of harmonic oscillators with 
the following Hamiltonian: 



B 



3 ^^^-^M 



2mj 2 



» ^ rrijUJjXjoXj 



^ V 2M,- 2 



(n2 - ni) + 



' ' ' ' {N1-N2). 



(10) 



Here {xj,pj} are positions and momenta of the oscillators coupled to the electron subsystem, 
whereas the variables {Xj,Pj} are related to the proton environment. The electron and 
proton surroundings are characterized by their own sets of effective masses rrij and Mj as 
well as by the two sets of eigenfrequencies Uj and Qj. The strengths of the couplings to the 
environments are determined by the shifts xjo and Xjq of the equilibrium positions of the 
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corresponding jth-oscillator. The bath Hamiltonian, Eq. (fTOl) . can be rewritten in the form 

v-/Pi , "^i^|[a;j + (l/2)x,o(n2-ni)]2\ 
^-^4^^ + 2 J + 

^(P! M,mX, + {l/2)X,,{Nr-N,)f 
^ \ 2M. 2 

\K{n^+n2)-\h{N^ + N2), (11) 

where the parameters Aa and are reorganization energies for the electron and proton 
environments, 

k^yP^- a.=5:Mi^. (12) 

3 j 

The systems of independent harmonic oscillators are conveniently characterized by the spec- 
tral functions Ja{uj) and Jbij^), defined as 

u^) = E -^-f^^i^ - ^.)' M^) = E -^-^^(^ - ^.)' (13) 

so that 

Xa= / —Ja{oo), Xb= / —Jb{oo). (14) 

Jo ^ Jo ^ 

D. Total Hamiltonian 

The total Hamiltonian of the system incorporates all the above-mentioned terms, as 
if = ifo + E ^kacl^Cka + E ^il^Kl^^il^ + ^Fai'a2&^ &1 + Vpataibtb2 - 

ka qP 
- E ^kRclnai - E ^InO'l'^kR - E ^kLclLa2 - E ^kL'A'^kL " 

J\y 

E ^PgNdg^bi - ^ T*^btdqN - E Tqpdgpb2 - E T*pbtdqP + 
9 g <? <? 

^ mja>|[xj + (l/2)xjo(n2 - 

^ ^ P,^ , M,l]f [X, + (l/2)X,o(iVi - iV^)]^ ^ 
^l2l^ + ^ J' 

where the Hamiltonian 

Ho = E(^-^- + ^-^-) - E 

u^„in„N„i + Uenin2 + UpNiN2 (16) 



is characterized by the renormahzed energy levels, 

e. = 6(°) - (1/4)A„ = - (1/4)A,. 

Here the repulsion potentials, Ue and Up, also incorporate shifts proportional to the cor- 
responding reorganization energies, Aa/2 and Af,/2. With the unitary transformation, 
tj = UaUb, where 

Ua = exp[-{t/2)Y,PjXjo{ni - na)], Ub = exp[-{t/2) - iVi)], 

3 j 

we can transform the Hamiltonian if, Eq. ( ITSl) . to the form 

ka qf3 
k k k k 

i:(|-'^)^i:(i^M). a^. 

where 

= (l/h) J2p, X,o, ^6 = (1/^) J] X,0, 

are stochastic phases operators, and ^ = + ^fe- The result of this transformation follows 
from the fact that, for an arbitrary function the operator U produces a shift of 

the oscillator's positions: 

U+^x,,X^]U = + {l/2)x,o{n, - n^), X, + (l/2)X,o(iV2 - N^)]. 

In addition, this transformation results in phase factors for electron and proton amplitudes: 

and 
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E. Combined electron-proton eigenstates and energy eigenvalues 

The electron-proton system with no leads can be characterized by 16 basis states of the 
Hamiltonian Hq: 

|1) = |Vac), |2) = a+|Vac), |3) = a^|Vac), |4) = 6+|Vac), |5) = 6^|Vac), 
|6) = a^6^|Vac), |7) = a+bt\Va.c), |8) = a+6+|Vac), |9) = a^6^|Vac), 
|10) = a+a^|Vac), |11) = a+a^6^|Vac), |12) = a+a^6^|Vac), |13) = 6^6^|Vac), 

|14) = a+6+6+|Vac), |15) = a^6^6^|Vac), |16) = a+a^6^6^|Vac). (18) 

Here, |Vac) represents the vacuum state, when both electron active sites and both proton 
sites are empty, whereas, for example, the state |7) = a^fej'lVac) corresponds to the case 
when one electron is located on the site le and one proton is located on the site 2p. The 
state |8) = 6^ I Vac) is related to the opposite situation with a single electron on the site 2e 
and one proton on the site Ip. It should be also noted that any arbitrary operator A of the 
electron-proton system can be represented as an expansion in terms of the basis Heisenberg 
matrices pj^ = \m){n\ {m^n = 1,..,16): A = J^mn-^rrmPm- We will also use notations 
Pm = Pm the diagonal operator. Thus, the operators {cti, ^2, ^2} can be represented 
as 

ai = pl + pl + pl + pl' + pl'+pl' + p\t + plt 

a2 = pI + pI + pI- pI' - pI' - pI' + p\l- p\l, 

bi=pt + pt + pl + P\l + Pf + P\' + Pf + P\l 

h. = p\ + pl + pl + p}^ - p'l - pI' - pf - pf,. (19) 

The Forster operator in the Hamiltonian if, Eq. (|T71) . given by aj'"a26^6i, is responsible 
for the electron transition from the electron site 2e to the site le accompanied by the 
simultaneous proton transfer from the proton site Ip to the site 2p. In the basis introduced 
above, the Forster process corresponds to the transition of the electron-proton system from 
the state |8) to the state |7) : ala2 6^61 = |7)(8| = p^. Using the eigenfunctions, Eq. ([T8l) . 
we can rewrite the Hamiltonian Hq in a simple diagonal form: 

16 

Ho = ^em pm, (20) 

m=l 
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with the following energy spectrum: 

El — 0, 82 — ei, £3 = 62, £4 = El, 

£5 — E2, — ei -\- El — uii, 
e-j = ei + E2- U12, Es = €2 + El - U21, 

£9 = 62 + -B2 - 1^22, £10 = ei + 62 + Ue, 
£11 = ei + £2 + -El - Uii - U21 + Ue, 

£12 = ei + 62 + -E2 - M12 - U22 + Me, 

£13 = £^1 + £^2 + Up, Su = ei + £1 + £2 - - U12 + Mp, 
£15 = e2 + -El + £2 - 1*21 - U22 + lip, 

£l6 = ei + €2 + -El + £2 — Mil — M12 — M2I — M22 + Me + Up. (21) 

For the Forster component of the Hamiltonian Hp, and for the Hamiltonian ifjir describing 
the direct tunneling between the sites Ig, 2e and Ip, 2p, we obtain the expressions 

Hp = Vf e'^ + y; e-^^ (22) 

and 

//dir = -A„ e-^" (p^ + + + p}^) - AHpI + p^ + p? + pit) e'^« - 

A,(p^ + P? + P^ + piD e^^^ - A* e-'^^ {pi + p^ + p^ + pi?)- (23) 

It should be noted that the operators Hp and ifair are non-diagonal. 



III. ELECTRON AND PROTON CURRENTS 

The transfer of electrons (protons) can be quantitatively characterized by the particle 
current flows between left/right (negative/positive) reservoirs, (//?), which are defined as 

k q 
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with indices a = L,R and (3 = N,P. Taking into account the equations for electron and 
protons amphtudes in the leads, 

i CkL = €-kL CkL - tkL a2 C^^", 

i CkR = ekR CkL - tkR €'2^" ai, 



i d, 



qN 



EqN dqN 



TgN hi 6 5 



i = Egp dgp " TgR 6 2 63 



(25) 



we obtain for the currents, 



iL = i^tkL{cl^a2e2^'') + h.c; iR = i^tkR{e 2«°c^^ai) + /i.e.; 

k k 

lN = tJ2 T<iN{d+^hie^^^) + h.c- Ip = iJ2 Tqp{e-^^^'d+ph2) + h.c. 
q <? 

It follows from Eq. ( l25l) that the leads' responses are described by the formulas 



(26) 



CkL = 4l - ^kL j dti gl^{t,ti) 02(^1) e 
dqN = 4/J - j dti gf^{t,ti) 61 (ti) e 



55a(il) 

16 (ii) 



5 



(27) 



etc., where 



9Lit, h) = -I e^--(*-*^) e{t - h), gj'Jt, ti) = e 



are the retarded Green functions of electrons and protons in the leads, c'kl-.d^^p are unper- 
turbed electron and proton operators in the electron reservoir a and in the proton lead /3, 
respectively, and 0{t) is the Heaviside step function. Within our model, we assume that 
electrons and protons in the leads are characterized by the Fermi distributions 



fai^^ka) 



exp ( 1 + 1 



1 -1 



, Eqf3 - fli3 , 

exp ( 1 + 1 



n -1 



respectively, having the same temperature T {ks = 1). However, the chemical potentials of 
electrons in the left (fiL) and in the right (/i^j) lead, as well as chemical potentials of the 
protons from the negative side of the membrane (/iA?) and from the positive one (/xp), can 
be different in the non-equilibrium case: 



^^L — ^^a + K, A^_R — ^^a■, /^a^ — ^^h■, /^p — /^t + K 



pi 
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where Ve and Vp are electron and proton voltage build-ups, and jj^b are equilibrium chem- 
ical potentials of the electron and proton reservoirs, respectively. Notice that the absolute 
value of the electron charge, |e|, is included into the definitions of voltages Ve, Vp, which 
are measured here in millielectronVolts (meV). Thus, the correlators of the unperturbed 
operators are given by 

= ^./^(^.-) e^^-(*-*^). (28) 

In the wide-band limit, it is convenient to introduce frequency-independent densities of 
electron (proton) states, 7q, (T/j), as 

7« = 27r ^ - eka); rp = 2nJ2 l^<?/3|'5(^ - Eg^)- (29) 

k q 

It should be noted that the currents ia and Ip are involved in the equations for the averaged 
populations derived from the Hamiltonian, Eq. f|T71) . 

{rii) = -iVF{afa2btbie'^) + iVp{e~'^ajaibib2) + iAl{afa2e'^'') - iAa(e"*^"a^ai) - in, 
(ri2) = iVF{ata2b^bie''^) - iVp {e~^^ aib^ b2) + zA„(e"*^"a^ai) - iA*(a|a2e''^°) - z^; 
(iVi) = zVF{ata2btbie'^) - z\/;(e-'«a+ai6+62) + lAKe-'H+bi) - iAb{b+bie'^') - In] 
(iVs) = -tVF{ata2btbie'^) + iV^{e-'^ataibtb2) + iAb{btbie'^') - tAl{e-'Htb2) - /p(30) 

Here, the brackets (..) denote averaging over the equilibrium states of electron and proton 
reservoirs, complemented by the averaging over fluctuations of both dissipative environ- 
ments. It is evident that in the steady-state regime, when the time derivatives of all popu- 
lations are zero, the electron and proton currents are determined by the Forster process and 
by the direct tunneling: 

= -i^ = iVF{ata2btbi e'^) - zV^(e"'^ ajaife^fea) + 
iAaie-'^"" a+ai) - iA*(a+a2 e^^"^), 
= -Ip = iVF{ata2btbi e*«) - iV;{e-'^ 0+016+62) + 

zAlie-'^^ 6+62) - ^A;,(6+6i e*^''). (31) 

We assume that the Forster energy Vp, the direct tunneling rates, A^ and A;,, as well as the 
rates 7q and F^, which describe the tunneling between the active sites and the reservoirs, 
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are small enough compared to a parameter v AT which defines a characteristic energy scale 
of the noise operator ^ = + ^b, with a combined reorganization energy 

X = Xa + Xb- 

Then, all calculations can be done with an accuracy up to second order in the Forster energy, 
|VfP, and up to second order for the direct tunneling rates, |Aap and lA^p. The electron 
(proton) current consists of two components, iaF (I/sf), related to the Forster process, and 
ia,diT (-^/3,dir), describing the contributions of direct tunneling to the electron (proton) flow. 
The Forster components of the electron and proton currents are given by the same expression 
(up to the total sign): 

iRF = -iLF = IpF = -Inf = tV^{e-'^pl) - iVf{p\ e'^). (32) 

The direct electron (proton) current iR^div (-^Af.dir) is proportional to the tunneling rate 
A, (A,) : 

lN,dn = -/p,dir = iAl{e-'^''{pl + pI + pI + p\l)) + h.C. (33) 
A. Calculation of the Forster current 

To calculate the Forster component of the current up to second order in the energy Vf, we 
derive the Heisenberg equation for the operator pj neglecting the coupling to the reservoirs 
and the direct tunneling: 

'j/7 = ^ P7 + yFe-'HP7-P8), (34) 
where 6 is the detuning between the electron and proton energy levels, 

S = es- ej = €2 - ei - E2 + El - U2i + mi2- (35) 
The solution of Eq. fl34l) . 

P?W = -^Vf f e-*^(*-*^)e-^«(*^)[p,(ti) - p,{h)l (36) 

J —CO 

should be substituted in Eq. (132!) for the current i^P) 

i^p = -\Vf\^ t dti e-*'^(*-*^)(e-*«(*^)e*«W)(p7 - p8){ti) + h.c. (37) 
J —00 
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Here, we separate the averaging of the environment phases ^ = ^a + ^b from the operators of 
the electron-proton subsystem. For independent electron and proton environments, when 

we can also calculate the electron and proton functionals separately. In particular, for the 
electronic environment characterized by the operator C,a = Ylj^joPj (from here on h = 1) 
we obtain the relation 

exp{-iUt)} expiiUh)} = exp{-i[Ut) - Uh)]} exp{{l/2)[Ut),Uti)]-}, 
where the commutator, 

(l/2)[^a(t),^a(ti)]_ = -i^ TTijUjjx'^jQ sin ujj{t - ti), 

j 

is determined using the free-evolving oscillator operators, 

Xj(t) = Xj(ti) cosujjit — ti) H — sinti;,-(t — ti), 

Pj(t) = Pjiti) cosu!j(t — ti) — nijUjXj sin Uj(t — ti). 

For the Gaussian statistics of the system of independent oscillators, the characteristic func- 
tional has the form 



with 



h[ut),uti)]+) = = 5^(p,')x2ocos^,(t-ti). 



'2 



Taking into account the expression for the equilibrium dispersion of the jth-oscillator mo- 
mentum, (pj) = (mjCt;j/2) coth(ct;j/2T), we obtain the well-known expression j23| for the 
functional (e~*^"*^*)e*^''(*i)): 

{exp{-iUt)} expitUh)}) = exp{-tWiait)} exp{-W2ait)} , (38) 

where 

^ia(^) = /, ^"^^smtUjt = / dco —smtot, (39) 

.2 Jq u! 
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and 



W2a{t) = 



coth 



1 — coso;,- 



•jt) = / duj — coth I — 1 (1 - coscjt). 

Jo oj \zi / 



(40) 



Similar relations between Wn,{t) ,W2b{t) and the spectral function Jhioj) take place for the 
proton dissipative environment. Notice that for this model, the effects of the electrons and 
protons on the environments are disregarded. In the semiclassical approximation (T ^ uj) 
and for slow enough fluctuations of the environments {ojt -C 1), the functions Wiait), W2a{i) 
have simple forms 

WUt) = Kt, W2a{t) = \aTt\ 



Thus, we have 



(exp{-i^a(t)} exp{z^a(ti)}) = exp{-zAa(t - ti)} exp{-AaT(t - tif]. 



(41) 



The total characteristic functional involved in Eq. (!37|) for the Forster current, 
^g-*?(t)g«?{ti)^ = Q-^Mt-ti) ^-\T(t-tif ^ Yi^^ effective correlation time {h = 1), 



1 



Tr. 



'XT' 



which is determined by the combined electron-proton reorganization energy, A = + A;,. 
At strong enough electron-proton couplings to the surroundings, the correlation time 
is much shorter than the time scale of the probabilities p„, so that in Eq. ( !37l) we can put 
{Pt — Ps)(ti) ~ {pf — ps) (t) . It allows us to obtain a simple expression for the Forster current: 



^RF — —'t'LF 



IPF 



(42) 



where k looks like the well-known semiclassical Marcus rate 

4AT 



23 



TT 



^I^.Pexp 



(43) 



but with the only difference that instead of the reaction free energy of a proton pumping 
step, AG ~ — -El ~ £2 — ci, here we have the electron-proton detuning, 



5 = e2 - ei - E2 + El - U21 + M12, 



which is much smaller and can be even zero for the case of an exact electron-proton resonance. 
Near these resonant conditions, when 6 = X, the proton pump should be most effective. 
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B. Direct currents 



Similar calculations (not shown here) demonstrate that the direct electron (proton) cur- 
rent, Eq. fl55]) . is proportional to the standard non-resonant Marcus rate ka {kb): 



-^Ar,dir = --?'p,dir = kb{p^ + P7 + P9 + Pl2 " P4 " P6 " PS - Pll), 



where 



Kb 



TT 



|A„|^exp 



le2 - ei 



TT 



I Afep exp 



4A,r 

{E2 — El — XbY 



(44) 



(45) 



The processes of direct electron and proton tunnelings lead to the downhill transfer of pro- 
tons, discharging the proton battery. However, this process is significantly suppressed when 
the separation of the proton energy levels is much higher than the reorganization energy A^. 



IV. DENSITY MATRIX 

The electron and proton currents, Eqs. fH21) and flUl) . are determined by the diagonal 
elements of the density matrix of the electron-proton system (pm) over the eigenstates, 
Eq. f|T8|) . of the Hamiltonian, Eq. f|T6|) . To obtain the diagonal elements of the density 
matrix, we write the Heisenberg equation for the operators pm taking into account the basis 
Hamiltonian Hq = 'Ylm^npni complemented by terms which are responsible for: (i) the 
Forster process Hp-, (ii) the direct tunneling events between the active sites ifdir, and (iii) 
the tunneling coupling between the reservoirs and the active sites -fftum 

ipm = [H,Pm]- = [pm, Hp]^ + [pm, Hdir]- + [Pm, Htun]- ■ 

With the tunneling Hamiltonian, Eq. ([3]), where the electron and proton operators are 
represented as expansions, 

(^cr / J Oio-;mn PrM ^o" / ^ ba;mn Pm 
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(see Eq. (fT9|) ), we obtain the contribution of the two pairs of reservoirs to the evolution of 
the operator pm as 

\pmi -f^tun]— ~ ^ ^ ^kR 6 ^ ^ ^kR ('^1;'"". Pm ^l;nm Pn ) 

^fci C^L (c^l;™?! Pm ~ 0'l\nm Pn) ^ ~ 
Y,T,N d+,N {KmnPl ' h.nm Pn) c'^^'^ - 
T,P e''^-'^ <p {b^,mn Pi - b2;nm Pn) ' {h-C}, (46) 

Substituting Eq. (1271) for the leads reactions, and averaging over the Fermi distributions of 
electrons and protons in the leads and over the fluctuations of the environments, we obtain 
the contribution of leads to the master equation for the probabilities (pm) '■ 

{[Pm, i^tun]-) = ^ 5^(7^ (Pn) - 7^ (P™)), (47) 

n 

with the relaxation matrix 

Imn = 7fi{|«l;mnni - /i?(^^nm)] + |cH;nmP/R(t^mn) } + 
7l{| I fL{uJm.n)} + 

^N{\bl-mn\'^[^ — -^Af(l^nm)] + | &l;nm P-^A^ ("^mn) } + 

^nm)\ ~r "2:nm 

YFp{Uran)}. (48) 

The products of free reservoir operators, such as c^'^(t), and an arbitrary Fermi operator of 
electrons, Zp, can be calculated using the formula 

(ZHt)cS(t)) = dt, {ctl\t,)ctl{t)) ([ZHt),a.(ti)]+)^(t-ti). (49) 

Similar formulas can be employed for the proton component. The Forster process contributes 
to the evolution of two components of the density matrix, pj and pg, 

[P7, Hp]. = -[p8. Hp]. = Vp p«7 - V* e-*« pi (50) 

Due to the weakness of the tunneling processes, we disregard the overlap of the different 
tunneling mechanisms in the master equation for the distribution (pm)- Substituting Eq. (15^ 
for the operator p^ and its conjugate jointly with Eq. fKTl) for the characteristic functional of 
the environments, we obtain the contribution of the Forster process to the master equation 
as 

{[p,,Hp]_) = -{WHp]_) = m((p8) - (pv)), (51) 
18 



where k is the resonant Marcus rate, Eq. (H3l) . In a similar way, we determine that the direct 
tunnehng between the active sites contributes to the equations for the following probabilities: 



(fp2, -f^riirl 


1 \L/ -J' uiij / " a\\t^o/ 


- (P2)), 


/r zj ^ 


-) = -([P5, ^dirj-) = ifib{{p5) - 


- (P4)), 


/ r Tj ^ 

UPe, -HdirJ- 


\ ■ { 1 \ 1 W \ ■ { 1 \ 

-) = if^aiiPs) - (Pe)) + if^biiPr) 


- (Pe)), 


[P7, -^dir]- 


_ = iHaiipg) - (Pt)) - il^biipl) - 


- (Pe)), 


[PS, -f^dir]- 


= -iKaiips) - {pe)) + inbiipg) 


- (P8)), 


[P9, -f^dir]- 


= -ifiai{p9) - {pi)) - if^biipg) 


- (Ps)), 


[Pll, -f^dir] 


- = -[pi2,-f^dir]- = iK,b{{Pl2) - 


(Pii)), 


[Pl4, -f^dir] 


- = -[Pl5,-f^dir]- = il^a{{Pl5) " 


(Pl4)), 



where ka and fc;, are the non-resonant Marcus rates given by Eq. Combining all con- 

tributions, we obtain the following master equation for the probabilities {pm)'- 

(Pm) + 7m (Pm) = X] 7mn (Pn) , (52) 
n 

with the relaxation rates ■jrn = '^n'^nm, where ■jrnn = Jmn given by Eq. for all matrix 
elements except 

72,3 = I2J + 73,2 = llj + ka, 74,5 = 7^5° + 75,4 = 7^ + h] 
76,7 = llj + A;;.; 77,6 = llj + h; 76,8 = 7m + 78,6 = 78?6" + ^a, 
77,8 = 77X + 78,7 = llj + 77,9 = IyJ + ka, 79,7 = 79 J + 

78,9 = 78!'9" + h; 79,8 = 79X + h; 711,12 = 7n!l2 + kb; 712,11 = 7i2"ii + ^'-s 

714,15 = 7uri5 + ka, 715,14 = 71m4 + ^a- (53) 

It should be noted that the key ingredient of the proposed model is the resonant Forster 
exchange of energy between electrons and protons. This process takes place in a time interval 

1 

'^F = — , 

2k 

where k is the resonant Marcus rate Eq. (H3|) . as follows from the solution of the rate 
equations, (py) = —n{p7 — Ps) = ^{Ps), derived in the absence of the leads. If our system is 
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initially in the state |8) with the excited electron and with the proton in the ground state, 
then, the probability to be in the state |7), where the proton is on the upper level and the 
electron in the ground state, is given by the formula 

Prit) = (1 - e~'-')/2. 

After a lapse of time scale Tp, the proton goes to the excited state with probability 1/2. 

V. RESULTS AND DISCUSSION 

The steady-state version of Eq. fl52l) . 

^7nm (Pm) = ^ 7mn (Pn), (54) 
n n 

{m,n = 1,..16), has been solved numerically jointly with the normalization condition 
TlmPm = 1) with subsequent calculations of the electron and proton currents through 
the system, Eqs. (H2l). fH^ . and populations of all active sites, {ria) and {No). To obtain 
numerical values, we assume that the electron potential well, presumably attached to the 
binuclear center, contains two active electron sites and has a radius tq of about 0.1 nm. The 
proton potential well with a radius Rq ~ 0.01 nm can be located at the pump center X at 
a distance i? ~ 1 nm from the electron sites. Thus, in a medium with a dielectric constant 
= 3 (dry protein), the Forster constant in Eq. ([7]) has a Vp ~ 1 meV. Taking into account 
renormalization effects for the direct Coulomb coupling between electrons and protons, we 
choose 

Mil — Ui2 — U21 — U22 = 400 meV 

which is close to the energy of the Coulomb interaction, u ~ 480 meV, of two charges 
located a distance i? ~ 1 nm apart. The on-site Coulomb repulsion energies,Me and Up, are 
estimated as 

Ue ^ Up 4000 meV, 

which is enough to avoid the double-occupation of the active sites. For the rates of the 
possible direct electron and proton transitions between the active sites, we take the values 
Aa = 1 meV and Af, = 0.1 meV, respectively. The tunneling couplings of the electrons to 
the leads are = Tpt = 0.85 meV, and the proton rates are F^r = Fp = 0.1 meV. For the 
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optimal efficiency of the pump, we choose the energy levels of the electron and proton active 
sites as 

ei = 100 meV, eg = 600 meV 

and 

El = 350 meV, E2 ^ 850 meV, 

ectron energy levels 62 and ei, corresponds to the realistic 



so that the difference between the e 
drop of the COX redox potential 
levels 

e2-ei = E2-Ei = 500 meV 



j2,[l^, and it is in resonance with the separation of proton 



We consider here intermediate values of the reorganization energies, 

Aa ~ Aft ~ 3 meV, A ~ 6 meV, 

which are higher than the Forster constant Vp and all other tunneling rates. Then the 
Marcus constants related to the direct tunneling, ka,kh, Eq. (H5l) . are negligibly small (~ 
IQ-wo nieV//i); however, the Forster rate, Eq. (1431) . is quite pronounced, k ~ 0.1 meV/h ~ 
150 ns^^. The rates Kb, and k can be measured in the units of meV/h or in the inverse 
nanoseconds (ns): 1 meV/h ~ 1500 ns~^. The real values of the reorganization energies 
Xa, Afe are not known yet for the enzyme cytochrome c oxidase, although it is expected 
that they are of order or higher than 100 meV 23 1. These numbers can be estimated 



from measurements of the temperature dependence of the Marcus rates K,a, f^b ()45|) for the 
transitions between the active electron and proton sites. 

It should be noted that at the reorganization energies A^, A^ ~ 100 meV, and at the 
physiological temperature, T = 36.6°C, direct tunneling processes are also significantly 
suppressed, 

Ka ~ 10"^ ns-\ Kb ~ 10"^^ ns-^ 

However, the Forster mechanism of energy transfer survives near the electron-proton reso- 
nance with the rate k ~ 30 ns~^. This means that even for the case of strong coupling to the 
dissipative environments, the pure electron-proton Forster exchange (with no leads) occurs 
over the time scale 

tf = 1/(2k) ~ 20 ps. 



21 



In the following, all contributions of the direct tunneling are disregarded, so that the 
total particle current is exclusively determined by the Forster component, Eq. fH2l) . and 
the electron flow from the left reservoir to the right one, in, is exactly equal to the particle 
current of protons, 

Ip = —In = ^R, 

flowing from the negative side to the positive side of the membrane against the concentration 
gradient. In other words, one proton is pumped through the membrane per each electron 
transferred to the oxygen molecule O2 that can play the role of our right electron reservoir, 
consistent with experimental observations of Refs. ^, 0, Ej. It should be mentioned that 
in the present model, we do not consider substrate protons, which are also taken from the 
negative side of the membrane to form the water molecules. 



A. Pumping effects 

Here, the positive direction of the current is defined to be from the higher chemical 
potential to the lower chemical potential. The electrochemical potential of the left electron 
lead, /iL, is chosen to be higher than the potential of the right lead at the positive voltage 

Ve-. 

HL = Ve, fiR = 0, 

whereas for the protons the chemical potential of the positive side of the membrane, hp, 
exceeds the potential of the negative side at the positive voltage Vp-. 

Hp = Vp, hn = 0. 

Notice that throughout the paper the "voltages" Ve, Vp incorporate the absolute value of the 
electron charge and are measured in meV. When the electron voltage is positive, Ve > 0, the 
electron particle current i^, Eq. (^^, should be positive because the electron concentration 
of the right lead increases. At normal conditions, the protons should also flow from the 
positive side of the membrane (having a higher chemical potential at > 0) to the negative 
side, so that the population of protons on the negative side should grow, that corresponds 
to a positive particle current In- 

In Fig. 2, we present the numerical solution for the dependence of the proton current 1^ 
on the electron (V^) and proton (Vp) voltages at the physiological temperature T = 36.6°C, 
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with E2 = 850 meV. The particle current is measured here in the inverse nanoseconds, ns~^, 
so that, for example, the value In = —1 ns^^ corresponds to the transfer of one proton per 
one nanosecond from the negative side of the membrane to the positive side. It is evident 
from Fig. 2 that the uphill proton current (corresponding to negative values of In) starts 
at electron voltages exceeding a threshold value Veo = 550 meV provided that the proton 
voltage build-up is less than 450 meV. At these voltages, the states 

|7) = a+6+|Vac) and |8) = a^6^|Vac) 

participating in the Forster transfer (see Eq. (142!) ) and having energies ~ 550 meV begin to 
be populated. It is of interest that at lower voltages the state |6) = a'^b'll'Vac) containing 
an electron in the state Ig with energy ei = 100 meV and a proton in the state Ip, having an 
energy Ei = 350 meV, is partially populated. Here, the electron-proton Coulomb attraction, 
Mil = —400 meV, comes into play, lowering the total energy to the value = 50 meV. 

For the chosen parameters, the particle current In saturates at electron voltages higher 
than 700 meV with the value corresponding to the translocation of 30 protons in one nanosec- 
ond. It shows the efficiency of the Forster pumping mechanism, although the real rate for 
the proton transfer through the D-pathway (see Ref. [3j) is much less: ~ 10^-10^ protons 
per second. This pumping rate can be obtained in the framework of our model if we sig- 
nificantly decrease the tunneling couplings between the active sites and the electron and 
proton reservoirs: ~ Tji ~ 10~^ meV, F^v ~ Fp ~ 10^® meV. It has no effect on the 
main features of the present model, and, in the following, we return to the case of the fast 
electron and proton delivery to the active sites. 

If the electron voltage is low enough, Ve < 300 meV, but the proton voltage is high, 
Vp > 500 meV, the proton flow reverses its direction, so that the protons move along the 
concentration gradient from the positive side of the membrane to the mitochondria interior. 
The downhill flow of the protons is especially signiflcant when the proton voltage exceeds the 
value of 850 meV. However, even at high proton voltages, the discharge of the mitochondrion 
battery can be prevented by applying the electron potential above the threshold VgQ = 550 
mV. We emphasize that, within this model, we do not need any additional gates to inhibit 
the translocation of protons back to the negatively- charged interior, although the pump can 
work in the reverse regime. The optimal value for the proton voltage build-up, Vp = 250 
meV, correlates well with experimental data for the proton-motive force of about 200-250 
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meV |2|,y,B. 

The resonant character of the Forster energy transfer is demonstrated in Fig. 3 where we 
plot a dependence of the proton current Jtv on the variation of the higher energy level of the 
protons, -E'2) at several temperatures T measured in degrees Celsius. It is evident that the 
current In has the maximum absolute value at the energy 

E2 = e2 - ti + El - \ = 844 meV, 

which is slightly shifted from its resonance value E2 = 850 meV in accordance with the 
maximum of the Marcus constant k, Eq. (H3l) . 

In Fig. 4 we present the temperature dependence of the uphill proton current near the 
optimal point 

Ve = 700 meV, Vp = 250 meV, E2 = 850 meV. 

It is clear that the proton pumping peaks at temperatures between 0°C and 100°C with a 
strong decrease when the environment is colder than the water freezing point 0°C. However, 
the effect survives much better at high temperatures. Curiously, for the parameters used 
the uphill proton current has a maximum at temperatures about that of the human body 
(36.6°C). 

VI. CONCLUSIONS 

In conclusion, we proposed and analyzed quantitatively a simple nano-electronic and 
nano-protonic model reflecting the main features of the electron-driven proton pump in the 
enzyme cytochrome c oxidase. We analyzed quantum-mechanical Hamiltonians for this sys- 
tem taking into account tunneling couplings of electrons and protons to their corresponding 
reservoirs and dissipative environments, as well as the electron-proton Coulomb interaction, 
including the resonant Forster term. Applying methods of condensed matter physics, we 
obtained expressions for the electron and proton currents as well as the equations of mo- 
tion for the density matrix of the system. These equations were solved numerically, and we 
demonstrated that the resonant Forster energy exchange between electrons and protons can 
lead to the proton transfer from the region with smaller proton concentration to the region 
with larger proton concentration, thereby achieving a proton pump. The dependence of this 
phenomenon on temperature and the system parameters were studied and we showed that 
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the proton pump works with maximum efficiency near physiological temperatures and at 
electron and proton voltage build-ups related to their values for living cells. 
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N: Negative side of membrane 




FIG. 1: (Color online) (a) Schematic diagram of the electron and proton pathways in cytochrome 
c oxidase with suggested locations for the active electron and proton sites, (b) Schematic energy 
diagram of the simultaneous electron and proton transport. 
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FIG. 2: (Color online) Proton current Ijv (a number of protons transferred through the membrane 
in one nanosecond) as a function of the electron (14) and proton (Vp) voltage build-ups at the 
physiological temperature T = 36.6°C and at the resonant condition, E2 = 850 meV. Notice that 
the absolute value of the electron charge, |e|, is included into the definitions of voltages Ve, Vp, 
which are measured here in meV. 
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FIG. 3: (Color online) Dependence of the proton current In on the resonant conditions (a variation 
of the upper proton energy level E2) at different temperatures, for optimal values of the electron 
and proton voltages: Ve = 700 meV, Vp = 250 meV. 
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FIG. 4: Proton current /^r as a function of temperature T for E2 = 850 meV, T4 = 700 meV, 
Vp = 250 meV. The maximum value of the uphill proton current |/Ar| (which appears as a minimum 
in the plot) corresponds to the temperature T = 36.6°C. 
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